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VORTEX INTERACTIONS AND DECAY IN AIRCRAFT WAKES 


Alan J. Bllanln, Milton E, Teske, Coleman duP. Donaldson, 

and Guy G. Williamson 

Aeronautical Research Associates of Princeton, Inc. 
Princeton, New Jersey 


SUMMARY 


The dynamic interactions of aircraft wake vortices are 
investigated using both inviscid and viscous models. The phenome- 
non of vortex merging which results in the rapid aging of a vortex 
wake is examined in detail. It is shown that the redistribution 
of vorticlty during merging results from both convective and diffu- 
sive mechanisms. 

Inviscid computations show that the merging phenomenon is 
sensitive to small changes in the spanwise load distribution. In 
particular, the fuselage vortex, shed from the wing-fuselage 
junction, can play a significant role in bringing about merging of 
wing tip and flap vortices. The usefulness of point vortex compu- 
tations in predicting merging is also explored, 

A vortex wake code including turbulent transport is developed 
using a second-order closure model of the turbulent transport. 
Vortex wake merging computations using this code indicate that a 
low hazard wake occurs when the generating aircraft trails flap and 
wing tip vortices from each side of the wing of the same strength 
and sign. To achieve this optimum, the flap vortex is located out- 
board at. approximately k0% of the distance to the tip vortex. 

Preliminary computations of multivortex wakes interacting 
with the ground suggest that the merging phenomenon may be altered 
by the presence of this ground plane. Computations have shown that 
jet exhaust does not alter the merging phenomenon. The effect of 
ambient atmospheric turbulence on the aging of an aircraft wake is 
investigated at constant turbulent dissipation rate. It is shown 
that under stable atmospheric conditions, when atmospheric macro- 
scales are less than or equal to the vortex spacing, misleading 
results may be obtained. 


1 . INTRODUCTION 


Low hazard aircraft vortex wakes are produced when the 
configuration of an aircraft Is such that multiple vortices are 
trailed which Interact strongly enough to result In merging. The 
mechanism of merging Involves the convective and turbulent redis- 
tribution of trailed vortlcity. The merging phenomenon In aircraft 
vortex wakes or the pairing , condensing, or coalescing phenomenon 
In turbulent shear layers is decidedly nonlinear. Thus, an 
experimental or numerical approach to the study of the phenomenon 
has been favored. There is now a need to understand it more fully. 
In particular, a model which could predict the redistribution of 
vortlcity during merging would aid us In our search for an aircraft 
configuration of low vortex hazard. The work reported here reviews 
the methods that are available to investigate aircraft vortex wakes 
and describes a new technique which has been developed to give a 
detailed description of the viscous interactions in a wake. This 
new technique makes use of a code which solves the fluid equations 
of motion Including turbulent transport through a second-order 
turbulent closure of the Reynolds stress equations. 

It is appropriate to introduce the merging concept here in 
order that the reader may Initially see the significance of the 
phenomenon. Consider the wakes of an aircraft in the cruise con- 
figuration and in the landing configuration as schematically 
illustrated In figure 1. The trailed vortlcity of the cruise wake 
is found in two regions and the distributions are nearly axisym- 
metric . Consequently, the production of turbulence and the subse- 
quent diffusion of vortlcity across the nearly axlsymmetrlc stream- 
lines are suppressed by the centrifugal effect of the swirling 
velocity-. The aging in the absence of atmospheric turbulence, wind 
shear, or stratification results only from the modest amounts of 
turbulence generated in the vortices. 

When an aircraft is in a landing or takeoff configuration, 
multiple vortices are trailed to form a wake. These vortices are 
both convected and strained by each other. This Induced straining 
field is responsible for the merging or pairing of the two vortices, 
wMch then results in the production of turbulence and even wider 
redistribution of the trailed vortlcity. We have tried schemati- 
cally to Illustrate in figure 1 the fact that at comparable times 
after aircraft passage the multiple-vortex wake has a more widely 
spread distribution of vortlcity than does the wake from an air- 
craft in the cruise configuration. It should be apparent that if 
circulation in the half-plane is taken as a measure of the strength 
of the wake, no reduction in the circulation can occur until the 
trailed vortlcity has been diffused to the wake centerline. The 
turbulence produced during merging provides the mechanism for this 
diffusion . 
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The viscous interactions in an aircraft wake cannot be 
ignored if there is to be any attempt to predict the decay of wake 
intensity with distance downstream of the aircraft. While this 
report treats both inviscid and viscous interactions, our emphasis 
is on transport phenomena. The report is organized in five parts; 
the first (Section 2) reviews inviscid techniques to compute wake 
structure and vortex-vortex interaction. The mechanism responsible 
for merging is discussed and an inviscid study of the vortex wake 
of a B-747 aircraft is described. In Section 3 a description of 
our turbulent model based on a second-order closure of the Reynolds 
stress equations is given. A description of the computer code 
which solves the turbulent model is given in Section 4, and a 
variety of vortex wake interaction problems are posed and solved in 
Section 5. Finally, in Section 6, conclusions and recommendations 
are offered. 
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NOMENCLATURE 

aspect ratio and turbulent model constant 
wing span and turbulent model constant 
tall span 
local chord 

sectional lift coefficient 

center of pressure 

wing reference chord 

rolling moment coefficient 

lift coefficient 

moment coefficient at = 0 

length 

gravity 

height of tall above initial wing wake vortex 
sheet; height above ground 

lift 

angular momentum in fluid in bounded region 
due to ith vortex 

pressure 

root mean square of twice the turbulent kinetic 
energy 

position vector 

vortlcity spread length scale 
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wing semi-span 
turbulent model constants 
reference area 
time 

velocity components in a Cartesian x,y,z system 

free stream velocity 

model diffusion coefficient 

Cartesian coordinates 

centroid of trailed vorticlty 

angle between local section of wing alid reference 
line when subscripted or standard angle of attack 
when not subscripted 

circulation 

Inboard trailing edge flap 
midspan trailing edge flap 
leading edge flap angle 
turbulent dissipation rate 

axial, streamwise, or x component of vortlcity 
fluctuating temperature 
mean temperature 

turbulent dissipation scale 
macroscale or Integral scale parameter 
kinematic viscosity 
density 

stream function 


follower or flap 

generator 

tall 

wing 
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2. INVISCID DESCRIPTION OP AIRCRAFT WAKES 


While a mathematical model of an aircraft wake must include 
turbulent transport if there is to be any attempt to predict aging 
of the wake, considerable insight into wake structure and dynamics 
can be obtained from an inviscid description of the wake. Inviscid 
techniques which are used to predict wake geometry and structure 
are reviewed here. Emphasis is placed on limitations of the models 
in an effort to make us cognizant of their range of applicability. 
The technique of using point vortices to approximate the rollup and 
subsequent vortex Interactions of an aircraft vortex wake receives 
special emphasis, and calculations are undertaken to systematically 
investigate the 8-7^7 at a lift coefficient of 1.2. What will 
become clear is that favorable vortex wake interactions are quite 
sensitive to small changes in wing load distribution. In parti- 
cular, inviscid analyses show that the strength and position of 
the fuselage vortex, which is trailed as a result of lift dropoff 
at the wing-fuselage junction, play a subtle roll in the favorable 
merging phenomenon. This result and other results which are 
obtained here bring out the need for measured wing load distribu- 
tion data. 


2.1. Method of Betz and Extensions 

The method of Betz (ref. 1) relates the load distribution on 
a simply loaded wing to two axlsymmetrlc distributions of swirling 
velocity in the distant rolled-up wake. The method is based on the 
application of certain '’exact" integral invariants of a two-dimen- 
sional distribution of vortlcity in a physically motivated manner. 
Many investigators have re-examined the Betz procedure and have 
extended the technique to include multiple vortices as well as 
distributed wing drag. These extensions have been recently 
reviewed in reference 2. 

The comparison of measured swirl velocity distributions in 
the wakes of clean aircraft (cruise configuration) with distribu- 
tions predicted by the Betz technique is remarkably good (ref. 3). 
However, when an aircraft is in a dirty (landing or takeoff) 
configuration, the Betz technique has predicted velocity distribu- 
tions which depart significantly from measured distributions (ref. 
4). This discrepancy has been explained by the observations first 
made by Dunham (ref. 5) of the strong Interaction and coalescing 
or merging of multiple vortices into a single, more diffuse vortex. 
The Betz method assumes a priori that each vortex trailed in the 
wake rolls up independently. That is to say that other vortices 
in the wake do not affect vortex roll-up or structure. The merging 
phenomenon is one in which the straining of vortex on vortex is the 
dominant interaction. Therefore, the prediction by the Betz tech- 
nique of the vortex velocity distributions in a multlple-vortex-palr 
wake, where vortices (particularly vortices of like sign, as we will 
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show) come in close proximity, must necessarily fall. The Betz 
technique has proven to be a simple means of obtaining a first look 
at the structure of aircraft vortex wakes. It cannot, however, be 
used in situations where straining interactions between vortices 
are important. The technique is still quite useful as a means of 
specifying approximate initial velocity conditions for the numer- 
ical wake code described in Section 4 . 


2.2 Invlscld Description of Merging 

A physical description of the merging phenomenon can be 
obtained from a simple model computation first undertaken by 
Corsiglla, et al. (ref. 6). If a weak vortex of vanishingly small 
strength is placed in the velocity field of a second vortex of 
strength V , the weak vortex is simply convected about the strong 
vortex. The structure of the strong vortex is unimportant s^nce 
it cannot change in this invlscld calculation. The weak vortex has 
a circular region of diameter d containing negligible vortlclty 
and is initially positioned at z/d = 1.0, as shown in figure 2. 

The strong vortex is positioned at y = z = 0 . As the computation 
proceeds, the weak distribution is strained about the strong vortex 
and, as shown, becomes progressively more and more wrapped about it. 
It is significant that the initially circular region becomes a 
spiralled sheet and that the fluid between the spiralled sheet is 
irrotational . Of course, transport will quickly act to smooth this 
vortlclty distribution so that the end state is one nearly axi- 
symmetric vortex with vortlclty distributed over a circular area of 
diameter greater than 2d . When the merging phenomenon occurs 
between vortices of nearly the same strength, the axial symmetries 
of both vortices are destroyed by the strain of vortex upon vortex. 
While the Initial phases of merging simply Involve convection, 
turbulent transport soon becomes the mechanism by which vortlclty 
is redistributed outward. When present in aircraft wakes, merging 
results in a significant reduction of the rolling moment Induced on 
an encountering aircraft. The outward redistribution of the trailed 
vortlclty, which promotes the rapid aging of the wake, is the 
mechanism responsible for the reduction in rolling moment. 

Since the initial phase of merging is a convective redistri- 
bution of vortlclty, attempts have been made to predict 
analytically and numerically the strain field which will not allow 
a vortex to remain an entity. Moore and Saffman (ref. 8) have 
given an analysis of a vortex of elliptic cross-section in which 
the vortlclty is uniform. This vortex is embedded in a flow field 
composed of a pure strain combined with a far field of solid body 
rotation; that is, the flow field Induced by one vortex on another. 
It has been possible to apply their analysis to determine the 
proximity d and strength f 2 at which a vortex of initial dia- 
meter d (3 and strength fi (uniform vortlclty) will not yield a 
steady vortex solution of elliptical cross-section. Unfortunately, 
this does not imply that merging is imminent, since unsteady 
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solutions may still exist. However, as we will show, numerical 
experiments using large numbers of point vortices to represent the 
distributions of vortlclty suggest that the criterion obtained is 
a good estimate of the strain required to bring about merging. 

The details of this computation may be found In Appendix A. 

In figure 3 are shown the curves of critical relative vortex 
strength versus vortex separation d/dc below which steady stable 
solutions do not exist. Also shown Is the merging boundary 
determined by Rossow (ref. 9) In his Investigation of the merging 
phenomenon between two vortices of uniform vortlclty having the 
same core diameter d^ . Typically, 4o point vortices are used to 
represent each vortex distribution. Also shown Is the critical 
separation ratio for r 2 /fi = 1.0 as computed by Roberts and 
Christiansen (ref. 10). Here It Is estimated that several hundred 
point vortices have been used to represent each distribution. 

The agreement Is reasonable. We point out, that in any case the 
analytic model Is not exactly the fluid flow being computed by the 
distributions of point vortices. That Is because the vortex In 
the analytic model which is providing the straining field Is not 
permitted to have any structure — it Is a point vortex of strength 
Tz ■ It Is, however, heartening to find the favorable comparison. 


Judging simply by the areas in which merging Is predicted, 
it Is suggested quite strongly that merging between vortices of 
like sign is more probable than merging between vortices of unlike 
sign. This observation has recently been confirmed In a wind 
tunnel using wings to generate tip vortices which are marked with 
smoke (ref . 7 ) • 


2.3. Wake Roll-Up and Vortex Interactions 
Using Point Vortices 

Point vortex computations have a widespread popularity among 
researchers and have been used to model shear layers as well as 
aircraft vortex wakes. The object of this section is twofold: 

a) to compute some measure of the roll torque experienced by 
a follower aircraft In the wake of a B-7^7 aircraft In trim config- 
uration modeled by a distribution of point vortices; 

b) to present a critical discussion in an attempt to define 
the practical limits of applicability of point vortex computations. 

Results obtained by the point vortex technique are compared with 
experimental data where possible. The procedure used to estimate 
the angular momentum which can result In a roll torque on a follow- 
ing aircraft that encounters the wake Is depicted In figure 4 and 
Is outlined below. 
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Figure 3. Merging boundaries for a vortex with a uniform vortlclty 
core as computed from the model given In reference 8 
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Basic Inputs are: 

1) Points specifying wing planform with 
no leading or trailing edge flap 
deflections 

2) Leading and trailing edge flap deflections 

3) Total lift coefficient 

4) Wing-tall areas and longitudinal positions 
relative to aircraft c.g. 

5) Desired number of chordwlse and spanwlse 
subdivisions of wing for vortex lattice 
model 

:: ~ ~ r 


Calculate spanwlse lift distribution, angle of 
attack, and tail lift coefficient. 


Desired maximum number 

"" 1 


of line vortices in 

— 



wake model 





Calculate strength and initial positions of 
wing and tall vortices. 


Calculate vortex trajectories as they convec- 
tlvely interact with one another. 


Using a hazard criterion such as angular 
momentum of the fluid in a bounded region, 
survey the wake and evaluate its hazard 
potential , 


Figure 4. Flow chart for wake hazard evaluation using discrete 
invlscid vortices 
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In Appendix B the procedures used to compute the 
discrete vortex wake behind a generating aircraft Is given. 

With the positions of discrete vortices given, the remaining 
step Is to evaluate how hazardous the flow field Is for smaller 
aircraft. ' - ■ 


Hazard criterion for vortex wakes modeled by distributions 
of point vortices .- Using strip theory, the roiling moment Induced 
by a point vortex located In the plane of a following wing of semi- 
span s Is proportional to 


<=£ dy (1) 

where the vortex Is located at y = d and the airfoil Is assumed 
to have constant chord and constant sectional lift coefficient. 

When d > s 3 the Integral Is well behaved; however, when -s < d < 
s , the Integral must be Interpreted as a Cauchy principal value. 

The difficulty, of course, occurs when d = ±s . Here the rolling 
moment Is undefined. This singularity has made it difficult in the 
past to interpret results of point vortex computations quantita- 
tively . 

The angular momentum In a circle of radius s Is a well 
defined quantity which Is related to the torque Induced on an 
encountering wing. Physically the magnitude of the angular 
momentum In the fluid Immediately surrounding the following air- 
craft must give a measure of the Induced rolling moment. Exception 
must be made for certain pathological situations which do not appear 
relevant to the discussion here. 

The angular momentum of the fluid In a circle of radius s 
due to a single point vortex as shown In the following sketch is a 
straightforward computation. The coordinate system is centered on 
the vortex center with y axis through the center of the circular 
window located at Tq . The angular momentum with respect to the 
center of the circular window due to the 1th vortex Is 
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( 2 ) 


= I (r “ r^) X V dA 
Substituting for each of the vectors 


/ 


? = yj + zk 


^0 = 


V = 


2TT(y^ + z^> . 


(zj - yk) 


and performing the Integration yields 


( 3 ) 

( 4 ) 

C5) 


^ (s^ - (6) 

It can be shown that vortices which are outside of the circu- 
lar window do not contribute to the angular momentum. The contri- 
bution due to each of the vortices In the wake model Is calculated 
and summed to find the total angular momentum of the fluid within 
the circle. The computation Is repeated with the circle moved 
laterally and vertically until the wake field at a particular down- 
stream station has been completely surveyed and the maximum value 
of angular momentum found. This value Is then used as a hazard 
measure for different wakes with different size follower aircraft. 

The techniques described above have been Implemented In a 
code and are now used to study the wake flow behind a B-747 air- 
craft . 


Results of Invlscld point vortex computation .- The Invlscld 
point vortex wake model was used to study theoretically the wake 
flow field generated by a B-747 aircraft at a lift coefficient of 
1.2. First the wing load distributions were computed for a matrix 
of flap settings. The c.g. position of the aircraft was assumed 
to be 29 % of the mean aerodynamic chord, and the aircraft was In 
trim with leading edge flaps fully extended — typical for landing. 
The following matrix of flap angle settings was Investigated and 
Is shown with the computed wing and tall coefficients. 


13 




A spanwise lift distribution for each corner of the matrix, 
as determined by the vortex lattice program, is shown in figure 5. 
At this point, we note what will be shown in the next subsection 
to be one limitation of the computed results presented here. The 
presence of the fuselage and the associated lift drop-off at the 
wing-fuselage junction cannot be calculated with this vortex 
lattice code. 

Each load distribution was discretized with approximately 40 
pairs of line vortices plus one pair for the tall vortex, and 
trajectories were computed to a downstream distance of about 
1.5^km. ^ Integration step size had to be a compromise between 
satisfying a Courant condition and avoiding accumulation of round- 
off error. In any event, a compromise was struck, and calculations 
all proceeded with the Integral invariants of Impulse and differ- 
ence kinetic energy (see ref. 2) not varying by more than a few 
percent from their initial values. 
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ure 5. Computed B-747 wing lift d 
extended 


VJl Pig 


istribution 



= 1.2, leading edge flaps 



Typical trajectories of the discrete vortices are presented 
and discussed In the next subsection, and we only present here the 
results obtained from our hazard criteria. 

Shown below Is the normalized maximum value of angular 
momentum in a circle whose diameter is -that of the Lear jet or DC-9 
aircraft, at a downstream distance of 1'.5 1cm. The normalizing 
angular momentum was the value computed with 6^.^ = - 30° 

(landing configuration) for the Lear jet follower. 



As a first general observation, the extreme variation of 
angular momentum over all configurations amounts to about 26 ^. 

This is disconcerting In that towing tank tests (ref. 5) have shown 
at these downstream distances rolling moment reductions In excess 
of 50^. This apparent discrepancy can be explained once we recall 
that point vortex computations, when carried out carefully, can at 
best only predict the convective redistribution of vortlclty since, 
by assumption, diffusion is neglected. This is not to say that by 
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using a suitably large number of point vortices the convective 
redistribution of mean velocity, as well as some turbulent struc- 
ture, could not be duplicated for some time. What is missing in 
the cascade energy transfer process is dissipation at the smallest 
eddy scales. The question which is difficult to answer is, 

"For just how long (time) can a point vortex computation be 
expected to yield a realistic simulation?" For surely, as time 
becomes large, the kinetic energy in the real wake cross-plane 
goes to zero, while a point vortex computation conserves kinetic 
energy of the fluid indefinitely. 

If we return to the question of numbers of point vortices 
needed, argument can be made that vortex placing should be of the 
order of the dissipation scale X which is a small dimension when 
compared to a relevant length scale in a vortex wake; say, the 
vortex spacing d . As a practical constraint, Rossow’s (ref. 9) 
results using 40 vortices per vortex in merging computations seem 
to compare favorably with our extension of Moore and Saffman's 
analysis (ref. 8) shown in figure 3. However, we note that Rossow 
has reported to us in a private communication that the reduction in 
wake strength after merging is not very significant. This- differs 
from the results of the full transport merging computations 
discussed in Section 

Returning to the normalized angular momentum matrix does show 
the general trend that proceeding from the standard landing config- 
uration, lessening the deflection of the midspan flap does, in 
general, reduce the angular momentum. The results for the 30°/0° 
configuration do not yield the lowest hazard wake for the Lear jet 
as has been observed in subscale and flight tests. However, the 
configurations bracketing the 30 °/ 0 ° configuration are the least 
hazardous. We speculate that the results for a 40°/0‘^ configuration 
would be even higher than the 30 °/ 0 ° result as more and more lift is 
carried on the inboard flap. Our calculations have shown that for 
the 30°/0° configuration the flap vortex is about 30 ^ larger than 
the tip vortex. The 20°/0° configuration has the flap vortex 
strength approximately equal to the tip vortex strength — a design 
criterion which is shown in Section 5 to be highly desirable in 
terms of maximizing vortex-vortex straining interactions. 

The absolute minimum hazard configuration for the Lear jet as 
identified from this computation is the 30°/10°, which is not 
easily explained at this time. We note only that in this case our 
computations have shown that the initial distribution of vortices 
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initially broke into three clusters, two of which then merged, 
leaving two distinct vortex pairs in the wake. The angular momen- 
tum contributions of the clusters were about equal in the 1.5 km. 
downstream position. 

Curiously, the most hazardous wake for the Lear Jet is the 
0°/30° configuration, a situation which has been observed by NASA. 

The results for the DC-9 (bf/bg = 0.5) encountering aircraft 
are, in general, not surprising. The angular momentum is, of 
course, greater since the DC-9 wingspan is greater than the Lear Jet 
wing span. Configurations which give hazard alleviation to the 
Lear Jet also, in general, are less hazardous to the DC-9. 


2.4 Selected Studies Using Point Vortices 

We have tried to stress in the previous section that point 
vortex computations do indeed contribute to our understanding of 
the complicated interactions occurring in vortex wakes. We have 
also, however, expressed caution in interpreting results, particu- 
larly in terms of assessing the deintensification of aircraft wakes 
as a consequence of merging. When done carefully, point vortex 
computations can at best only model the convective redistribution 
of trailed vorticity. By assumption, diffusion of vorticity and 
dissipation of kinetic energy are neglected. Nevertheless, point 
vortex computations are valid over time scales in which diffusion 
can be neglected. In terms of aircraft wakes, a practical time 
limit is up to the time of vortex merging. 

In this section, point vortex computations are used to obtain 
results which are relevant to vortex merging and wake interactions 
with the ground. 

If given that a vortex wake in the immediate vicinity of a 
generating aircraft is composed of a tip and flap vortex, the like- 
lihood of merging should Increase if the separation between flap 
and tip vortices does not Increase with distance downstream. 

Figures 6 and 7 . illustrate the trajectories of two differently 
structured wakes. The tight orbiting of flap and tip vortices in 
figure 6 will result in merging, while the weak Interaction 
between flap and tip vortices shown in figure 7 may not result in 
merging. 

For two-vortex-palr wakes, a qualitative prediction of 
whether merging may in fact occur can be obtained from a wake 
classification chart shown in figure 8 . When vortices are of like 
sign and fall in the region marked "remain together," one might 
expect that merging will eventual ly occur. On the other hand, 
wakes which fall in the region denoted "separate" may possibly 
remain unmerged. Wakes from real aircraft in a landing or takeoff 
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Figure 6 Vortex centroid locations as seen from down- 
stream with strong Interaction between 
neighboring vortices (ref. 2). 
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Figure 7. Vortex centroid locations as seen from down- 
stream with weak interaction between like- 
signed vortices -pairs diverge. The time 
interval between consecutive Integers is a 
constant (ref. 2 ) . 
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Figure 8. V/ake classification chart for two-vortex-pair wakes (ref. 2) 



configuration, particularly with aircraft having a split flap 
arrangement, have far more structure than can be represented simply 
by two vortex pairs. In fact, the B-7^7 aircraft in a landing 
configuration trails six discrete pairs If the tail vortex Is 
Included (ref. 2). 

At the next level of wake complexity, we consider a third 
vortex pair which can be used to model a fuselage vortex. We 
choose as a representative case the wake from the Landlng/0 config- 
uration studied by NASA. The load distribution shown In figure 9 
has been taken from the work of Clffone and Lonzo (ref. 12). Their 
code does attempt to predict loss of lift, at the wing- fuselage 
junction. The Betz rollup recipe will predict rollup Into three 
discrete vortices (ref. 2). The trajectories as viewed in a cross- 
plane are shown in figure 10. As can be seen, the distance between 
the flap and tip vortices continually increases downstream. 

Clffone and Lonzo conducted flow visualization experiments in a 
towing tank on this configuration and reported that merging did not 
occur. It is instructive to compute the separation between tip and 
flap vortices as a function of downstream distance while varying 
the strength of the fuselage vortex (fig. 11). By weakening the 
fuselage vortex, it is possible to reduce the separation between 
flap and tip vortices. The minimum separation between flap and tip 
vortices occurs when Fj - -0.3Ff . Further reduction in fuselage 
vortex strength serves only to again Increase the distance between 
flap and tip vortices. 

In figure 12, the sensitivity of distance between flap and 
tip vortices is Investigated as a function of position of the fuse- 
lage vortex. By moving the fuselage vortex from y^/s = 0.125 to 
0.2 (a shift of 7-5^ of the semispan), the flap and tip vortices 
can be brought to a distance of only 9^ of the wing semispan. 

The above computations serve to demonstrate the sensitivity 
of wake geometry to modest changes in wing load distribution. It 
is an unfortunate fact that techniques to give accurate predictions 
of load distributions are somewhat lacking, particularly for wings 
with flaps and slats deployed. The magnitude of the loss of lift 
at the wing-fuselage junction quite strongly influences wake 
geometry. There is now a need to know the wing load distribution 
more accurately than ever before. 

The problem of multiple vortex interaction in the vicinity of 
the ground plane is examined next. It is during landing, inbound 
from the middle marker to touchdown, that is the most critical time 
in which to avoid a hazardous vortex encounter. Any vortex allevi- 
ation scheme will, therefore, have to operate in ground effect. In 
figure 13 are shown the computed trajectories of vortices in the 
wake of an aircraft having split flaps with only the inboard flap 
deployed. Vortices trailed from each half of the wing interact 
strongly and remain in very close proximity. In fact, the tip and 
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Figure 13. Vortex trajectories from a LDG/0 configuration 




flap vortices which are Initially separated by 0.55s move to 
within about 0.08s downstream. If the fuselage vortex Is simply 
removed from the calculation, as shown In figure l4 , the tightly 
wound vortex pattern Is altered so that the distance between tip 
and flap vortices Increases downstream. This again warns that the 
fuselage vortex can play an Important role In the merging process. 
Now If the aircraft Is moved to an elevation of 0.5s above the 
ground and the fuselage vortex replaced, the vortex trajectories 
shown In figure 15 result. The results of these computations are 
summarized In figure l6. The above computations should warn 
that alleviation schemes Involving vortex merging should be tested 
in ground effect to make certain that ground-plane interactions 
are not adverse to the merging phenomenon. 

Before closing this section, we carry out several computations 
which suggest when It Is appropriate to replace a distribution of 
point vortices modeling a discrete distribution of vortlclty with 
one point vortex. 

Figure 17 shows the modified Landing/0 configuration 
discussed in ref. 12 at a Cl - 0.8. This load distribution was 
represented by forty discrete vortices of equal strength. The sub- 
sequent positions are shown in figure l8. Plotted also are the 
positions of the flap, tip, and fuselage vortices as obtained from 
a three-vortex-palr representation. These clusters are in reason- 
able agreement with the three-vortex-pair computations. 

For a second computation, we remove the fuselage vortex from 
the modified Landlng/0 configuration as shown by the dashed line 
on figure 17 j still using forty equal-strength vortices to repre- 
sent the trailed vortlclty. In figure 19 ^ a comparison of the many- 
vortex representation Is made with a two-vortex-pair model. The 
results are quite acceptable. 

Figure 20 shows the distance between the flap and tip vortex 
centroids, as computed above, with quite good agreement. The study 
shown in figure 12 Indicates that moving the fuselage vortex to 20^ 
of the semispan of the wing causes the flap and tip vortices to 
move into close proximity. The load distribution for the LDG/0 
configuration (fig. 11) was altered Inboard to create a fuselage 
vortex of strength Fj = -0.47Ff_. As It turned out, the fuselage 
vortex centroid was located at y/s = O .185 . The separation 
between tip and flap vortices as a function of distance downstream 
is shown In figure 21. For this configuration, minimum separation 
of d/s = 0.175 Is reasonably well predicted. However, the 
distances downstream at which this minimum Is achieved do not 
compare favorably. When the fuselage vortex Is removed so that 
flap and tip vortices move apart, the computation using a distribu- 
tion of vortices again compares favorably with computations using 
one point vortex to represent each discrete vortex in the wake. 

The discrepancy when distributions move together lies in the Incor- 
rect approximation that the centroid of a given distribution of 
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Figure l8. Comparison between using forty point vortices to 

represent the trailed sheet of vortlclty and using 
three point vortices located at the centroid of the 
tip, flap, and fuselage vortices as shown In 
figure 17. ^ denotes centroid computed from the 

three-vortex computation. 
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Figure 19, Comparison between using forty point vortices to 

represent the trailed sheet of vorticity and using 
two point vortices located at the centroid of the 
tn.p and flap vortices as shown in figure 17 (the 
fuselage vortex has been removed)* S c 3 enotes 
centroid computed fr*orn the two-vortox computation. 
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vortlcity moves at the velocity Induced at its centroid by other 
concentrations. In reference 2 we show that the velocity of the 
centroid of a bounded region of vorticity in an infinite fluid is 
given by 


d| 

dt “ r 


(7) 


where Uq is the velocity field induced by the presence of all 
other concentrations of vortlcity in the fluid other than the one 
whose motion is being computed. As two distributions come close 
together, Uq may vary significantly over ^ , and the motion of 
the centroid is not simply the fluid velocity at the centroid 
Induced by all other concentrations of vorticity. 

It has been our purpose in this section to point out that 
simple discrete vortex computations can be used to assess quickly 
whether particular lift distributions will result in a multiple 
vortex wake which might merge. When vortices of like sign remain 
in proximity to each other, the merging process must eventually 
occur. An estimation of the deintensification as a result of 
merging requires the use of a transportive fluid model, such as 
that presented in Section 3. In sum, Invlscld computations indi- 
cate that small changes in aircraft configuration and, hence, wing 
lift distribution can have a profound effect on wake structure. 
More accurate knowledge of lift distributions on commercial Jet- 
liners, in the landing and takeoff configurations, with and with- 
out landing gear deployed, would aid greatly in understanding the 
dynamics of aircraft vortex wakes. 


3. THE TRANSPORT MODEL 


A second-order closure of the Reynolds stress equations has 
been under development at A.R.A.P. for many years. The current 
modeled system of equations Is one which has been used to predict 
successfully both quantitative and qualitative features of complex 
turbulent flow fields (refs. 14-17). In this section, we review 
the Invariant technique of turbulent modeling. The reader can 
find a more comprehensive development in references l6 and 17. 
Vortex wake calculations solving the equations developed here will 
be presented in Section 5 of this report. 


The Navler-Stokes equations govern the solution to well-posed 
turbulent flow problems which are of interest In vortex wakes. 
Unfortunately 5 full numerical unsteady solutions of these equations 
are beyond the current technical computer capacities. Therefore, 
we try to obtain solutions for the ensemble-averaged mean field and 
ensemble-averaged second-order correlations. 

We begin by writing every flow variable with a mean and 
fluctuating part (e.g., ]J = + Ui). We then make formal substi- 

tution into the equations of motion and ensemble average to obtain 
the mean flow equations . 
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We have assumed that the flow Is Incompressible but have 
permitted small changes In density due to changes in temperature. 
The only effect of density variation Included is In the gravita- 
tional body force term in the momentum equation (ref. l6)(the 
Boussinesq approximation (ref. l8)). In vortex wakes, this body 
force term adds buoyancy for hot engine exhausts. The problem 
with eqs . (8)-(10) is that our decomposition into mean and 
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fluctuating parts has led to the addition of unknown correlations 
UiUj and Uj_0 which makes the system of equations incomplete. 
The'^ essence of first-order turbulent closure is to model these new 
terms as functions of first derivatives, e.g.. 


3U. 

I3f7 


( 11 ) 


Where K is an appropriate "eddy viscosity coefficient." This 
technique closes the system of equations but, in a great many cases, 
places a large, uncompromising burden upon the selection, use, and 
form of K — certainly rarely a constant. 


Recognizing this probl em, we endeavor to find "exact" 
expressions for and U 40 . This can be done by subtracting 

the mean equations and (10) from the complete equations, leaving 

equations for Uj_ and 6 which can be premultiplied and averaged 
to yield the proper exact expressions (ref. 16). 
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Unfortunately, this set still does not permit closure because of 
the presence of third-order velocity correlations •* 
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pressure correlations: 
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and dissipation terms: 
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In fact, for every additional new equation written, we would intro- 
duce successively higher-order turbulent correlations. Second- 
order closure stops here and attempts to model the above-noted 
expressions in some rational manner. That many other closure model 
schemes exist goes without saying (refs. 19, 20). 

We strive for an ’’invariant" model. Our definition may be 
taken in two ways: one, any expression that replaces any of the 

above-noted expressions must exhibit the same tensor symmetry and 
dimensionality as the replaced term; two, our final simple model 
with its inevitable modeling constants must be consistently and 
fairly applied without modifications to a wide variety of flow 
geometries and initial conditions. 

We choose to replace these higher-order correlations by 
appropriate second-order correlations and their gradients, mean 
flow gradients, and proper scaling and modeling constants. The 
velocity correlations, eq. (15) can be diffusive in nature and 
are therefore replaced by gradients of the a ppropri ate second-order 
correlation. Thus, for the term involving ^k^l^j » write 
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whei’e we have Introduced the length scale A and velocity scale 
q to make the substitution dimensionally sound. The scale A 
(the macroscale) measures the breadth of the turbulence under study; 
its description would generally vary from flow application to 
application. The velocity scale q is taken as the root-mean- 
square of twice the turbulent kinetic energy 

q = u^u^ = uu + vv + ww (19) 


The constant v^ adjusts any magnitude difference between the 
modeled and unmodeled terms and is the "diffusion coefficient" of 
the Invariant model. 

The pressure correlations, eq. Cl6) may be rewritten to give 
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The first two terms on the right-hand side are diffusive and are 
absorbed into the velocity correlation model, eq. (l8). The 
remaining term is a tendency-towards-isotropy term whereby the 
fluctuations attempt to redistribute turbulent energy between the 
three energy components. Following Rotta (ref. 21) 
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The dissipation terms are modeled so as to be isotropic, so that 
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where b is an additional model constant. 

Upon substitution of models, eqs . (18)-(21), the complete set 
of Reynolds stress equations becomes 
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Further considerations have prompted us to add the two extra model 
constants A (in the tendency-towards-lsotropy term of the velocity- 
temperature correlation equations) and s (in the dissipation 
terms of the temperature-temperature correlations). The fundamen- 
tal velocity scale q may be determined by summing eq. (22) for 
1 = j = 1,2,3 to obtain a differential equation for q^. The 
fundamental length scale A may be specified by some algebraic 
relationship peculiar to the flow under consideration, but gener- 
ality comes by writing a differen tial sc ale equation and actually 
solving for A as we solve for u^uj , Uj^0 , and 0^ . The 
process is admittedly more difficult because the scale equation 
does not follow directly from the equations of motion. However, 
assembling the terms most normally felt to Influence the scale 
results in the equation used at A.R.A.P. (ref. 17): 
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The first term on the right-hand side is the turbulent production 
of scale; the second term results in dissipation; the third and 
fourth terms are diffusion; and the fifth term is a production term 
which results from buoyancy. 

The modeled terms in eqs. (22) -(25) are not expected to 
contain all possible available Information present in the exact 
third-order correlations they replace. It would be Incorrect to 
believe that the Invariant model is more than a careful attempt at 
modeling the essence of turbulent fluid dynamics. We do not claim 
to solve the turbulence "exactly," only "accurately." The distinc- 
tion is very real; the value of any model rests with the ability to 
give a good description of the turbulent process in a wide variety 
of applications. That is where the modeling constants come in and 
where the second half of our "Invariant" definition is applied. 

With a model of this complexity, surely any single flow can be used 
to determine our constants (disregarding the equations for A ) 
such that that one flow is precisely described. However, the 
chances are high that those same constant values will not be very 
accurate for a different flow application. But it is our Intent to 
match a number of rather simple flows "well" and thus yield a fixed 
(or "invariant") set of values for the unknown constants A , b , 
s , and Vq that will match all of these flows respectably. With 
these constants determined, we can then make predictions in more 
complicated flows, knowing that our constants will permit a reason- 
ably accurate result. Our success to date verifies this approach 
(refs. 14, 22-24 ) . 

The most sensitive model constant is b . In a sense, this 
is good since we must pinpoint the value of b closely or else 
little will match anything. Analysis of homogeneous turbulence 
gives values of b between 0.08 and 0.12. Comparing with other 
data as well, Donaldson (ref. l6) chose 

b = 0.125 

Computer studies of simple axisymmetrlc jets and wakes (ref. 24) 
yield 

= 0.3 

The constant A can be shown to be equivalent to the turbulent 
Prandtl number and, for neutral atmospheric shear layers (ref. 25), 
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its value Is approximately 

A = 0.75 

By comparing the frequency spectrum of the temperature fluctuations 
with the spectrum of the turbulent kinetic energy , we find (ref. 

22 ) 


s = 1.8 

The values of the constants of the scale equation are determined by 
similar examination of the simple flows of axisymmetrlc jets and 
wakes, boundary layers, surface layers, grid turbulence, etc., to 
give (ref. 23 ) 

= -0.35 
s^ = -0.6 
S3 = 0.375 

Sj| = 0.8 

These constant values have been used to predict many simple flows. 
Including flows not used to set the constants. They have demon- 
strated the versatility of the invariant scheme and have consist- 
ently given us cause for optimism whenever a new flow situation or 
experimental result permits further verification. 

Before proceeding with our model results for vortices, we 
will outline the numerical procedure used to solve the partial 
differential equations presented here. Then in Section 5 we 
discuss the results of our vortex study to date. 
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DESCRIPTION OF THE PROGRAM 


The partial differential equations of motion developed in the 
last section are solved numerically to obtain solutions for vortex 
interactions and merging. -An extensive three-dimensional code was 
developed through modification and expansion of a basic code 
previously developed under contract to the Navy for studying sub- 
marine wakes (ref. 14). A fairly complete discussion of the numer- 
ical technique involved and the operational aspects of the program 
is contained in references 26 and 27. Here we wish to summarize 
the techniques employed in going from the differential equations to 
the actual solutions they produce. 

We first transform the partial differential equations into 
finite difference form. This step requires decisions regarding the 
computational grid or mesh and the form of the first and second 
derivative operators appearing in many places in the complete set 
of equations. Our code has been written so that it is forward- 
differenced in the marching direction, so that for either t or 
X , the marching derivative is written typically as 
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3x 


W - W , , 
new old 
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It is also center-differenced on a uniformly spaced grid in both 
cross-plane directions, so that for a typical derivative we write 

W - W 

9W y+Ay y-Ay 

9y " 2Ay 

where Ay is the mesh spacing and value of W at 

position y ± Ay . Second derivative differences are, in fact, 
derivatives of the first derivative difference. The general finite 
difference formulation on unequally spaced grids has been presented 
elsewhere (ref. 28). We merely reduce its complexity by working 
with uniform meshes. What this reduction does is give us a 

system that is second-order accurate in the cross plane. This 
feature yields dividends in solution accuracy. We use a direct 
solver for the solution of the pressure field (ref. 29) and this 
particular code requires uniform spacing. 

If the equations appear complex, in their differential form, 
they are much more complicated in finite difference form. Never- 
theless, an examination of the equations and our differencing 
scheme shows that at any point in the computational domain the 
unknowns there will be related to each other by the coupling of the 
variables through the production terms and dependent upon values of 
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the variable at all four perpendicular points Immediately adjacent 
to the point under observation. Several standard solution techni- 
ques are available for solving these equations. Certainly an 
explicit scheme exists where we march In x (or t) by evaluating 
all of the terms In the equations with the exception of the marching 
value (e.g., ^new^ with current, old (and, hence, known) values. 

The difficulty here lies In a restrictive step size Ax = Xnew"^old 
that must be imposed to keep the solution scheme stable. 


We choose, rather, to use the implicit algorithm applied to 
two-dimensional cross-plane flows (ref. 30). The scheme works by 
first linearizing the finite difference equations (of necessity 
using old values In some places) and then solving for, typically, 
Wnew t>y assuming that wherever W appears In its individual 
equation we write ^new (^^^^^^ber, in an explicit scheme, wherever 
W appears, we almost always write Woi^). We then apply a fairly 
standard algorithm to solve for Wnew all of the other 

variables at • This alternating direction, implicit, trl- 

dlagonal algorithm requires the storage of intermediate values; we 
sweep twice in each direction for each marching step size Ax to 
arrive at Xj^ew our solution there. 


The linearization is an outgrowth of the refinement of our 
numerical skills In solving turbulent flows. We break out every 
equation individually so that the only unknown (and, hence, 
implicit) variable possible in each equation is the variable for 
which the equation Is solved. Thus, although the W equation 
contains U , V , and P , we substitute 3 ^old a ^old 

wherever appropriate. In the convection term, we write 


W^ = W ^^W 

old new 


to keep the system linear. We also write the two shear-like terms 
as 


vv 


8W 

""y 9y 


+ c 


vv 


vw 


-e 


z 


9z 


+ C 

vw 


with Gy , Eg; , Cyy , and Cy^ evaluated at ( )old * rest 

of the '^W’s that appear are written at W^qw and, hence, unknown 
before application of the algorithm 

The set of finite difference equations is solved in a compu- 
ter program called WAKEX operational at the NASA Langley CDC 6600 
computer facility. The User Manual pertaining to program operation 
is in reference 26, Here we summarize the technique used to ini- 
tialize a computation; the next section of this report gives results 
for a variety of vortex flows which have been studied. 
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The code has been developed to simulate the merging and inter- 
action of a trailed vortlcity from a wing. With a given load 
distribution, a vortex sheet is trailed from the trailing edge of 
the wing. In principle, the velocity field associated with this 
trailed sheet could be input as an initial condition to WAKEX. 
Unfortunately, an accurate numerical resolution of this flow field 
would require a computational grid with a prohibitively large 
number of grid points. Rather, we input a Gaussian of vortlcity 
to represent each vortex in the wake. We have chosen to do this 
rather than use the velocity distributions given by the Betz model 
for two reasons. The first is that the Betz swirl velocity profiles 
are Invlscld and, therefore, do not have a viscous core. If these 
profiles were input as initial conditions, many grid points would 
have to be used to resolve the high swirling velocities. The second 
reason is that much computer time would be used simply to adjust the 
invlscld profiles to their viscous counterparts. The Gaussian 
distribution of vortlcity we use is chosen to have the same circula- 
tion centroid and dispersion as the vortex it replaces. 


As an example, the initial conditions which would be used to 
compute the decay of a linearly loaded wing are now given. If the 
load distribution on a wing of semi-span s is given by 
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( 26 ) 


and, for the linearly loaded wing, is r^ ^ 0.3s . 
vortex of strength Fq is positioned at y = 0.5s 
velocity distribution given by 


A positive 
with a swirling 


V = 



(27) 


where r^ = (y-y)^ + (z -z)^ . The image vortex at y = -0,5s 
is automatically accounted for through reflection boundary condi- 
tions at y = 0 . 


For an "interior” or flap vortex, the dispersion formula is 
modified as shown in reference 2 to read 
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and the interior vortex is positioned at its appropriate centroid. 

Because of computer core size and execution time costs, we 
have developed the program to permit no more than four such 
Gausslans in the computational grid, which itself can be no larger 
than 65 ^ 65 points. When an image plane exists (say, the other 
half of the wing), the program is organized to solve for the right- 
hand half of the flow and provide appropriate matching/reflection 
conditions at y = 0 to handle the image vortices present in the 
left-hand plane. The Initialization procedure Involves telling 
WAKEX where at (y^,Z(n) the vortices of strength and spread 

Vq are located. We have found that typically 5 grid points are 
needed in the e-foldlng of the Gaussian to give reasonable results. 
This would fix the computational grid spacing at about 0.05s for 
the above-noted example. 

It makes little sense to extend the computational domain far 
beyond the region of interest, yet we cannot set V = W = 0 along 
such boundaries, since every vortex pair has a dipole far-field 
and every isolated vortex has a swirling velocity which decays as 
r"^ . We know, however, that well outside the region of vorticlty, 
the vortex-induced velocities should behave, to lowest order, as 
V a z/r^ and W y/r^ . These expressions are appropriate for a 
single vortex but are also just the first terms in a far-field 
expansion of the moments of vorticlty distribution in the computa- 
tional domain. We have carried out the moment expansion to order 
2 where 

Inm " // dy dz ; n,m < 2 

and have included the appropriate Inm factors in the far-field 
boundary conditions on V and W . This analysis is summarized 
in Appendix D. By invoking this option in the program, we permit 
any accurate vortex representation to serve as velocity initial 
conditions in WAKEX. 

In the spirit of specifying equivalent Gausslans of vorticlty, 
we permit Initialization of temperature (or concentration) 0 , 
axial velocity U , and turbulent kinetic energy q^ by 

2 2 

Value = Maximum exp(-r /r^) 
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where r and r carry the same definitions as before. Thus, 
for any iet (of up to four each ) of 0 , U , or, , a specifica- 
tion of maximum value, spread rate, and location (y ,z ) of the 
center will permit Initialization of these quantities. 

Having specified the initial conditions on the appropriate 
dependent variables, the initialization is complete except for 
program execution configuration and printout specifications. These 
are covered in the User Manual. Perhaps, before summarizing the 
results of our study, it would be best to caution the reader 
regarding program performance. Even though we have had a surpris- 
ingly high success rate, that rate cannot be duplicated by others 
without following some rather basic rules. Perhaps all of the 
caution can be contained in the warning that ^WAKEX will not 
return an accurate solution except when the initialization proced- 
ure gives an accurate, compatible, and acceptable set of profiles. 
Enough points must be present in the computational grid to permit 
the complete equations to keep an accurate track of solution 
behavior. The boundaries of the computational domain cannot be so 
close to the region of variability (the yc>Z 0 values) that they 
give a rapid Jump to the far-field conditions. Nor can the solu- 
tion scheme be driven too fast in its simulation. We are generally 
restricted to permitting about 5 ^ change in the variables in the 
marching direction, while cross-plane velocities and step sizes 
must be compatible with the Courant condition 
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which does not permit the solution to march faster than it takes 
Information to convect from one mesh point to its neighbor. The 
initial stages of vortex roll-up and Interaction are convective in 
nature and demand that these criteria be followed. Once turbulent 
diffusion begins dominating the flow, the nearness of the boundar- 
ies will generally decide the overall time of the computer 
simulation. 

In the next section of this report, we use the Vortex Wake 
program to follow the behavior of various configurations with and 
without image planes for vortices, turbulent energy, axial velocity, 
and concentration initialized in ways to give an overall view of the 
phenomena of vortex merging and interaction. Our examples do not 
cover the full breadth of cases possible (nor, in fact, cases we 
have run) but are used to illustrate the power of the turbulent code 
in simulating vortex dynamics and to suggest the minimum hazard 
vortex wake configuration. 
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5. COMPUTATIONS USING THE "WAKE VORTEX" TRANSPORT CODE 


In this section, a variety of computations are presented to 
demonstrate the flexibility of the Wake Vortex transport code. 
Results of numerical merging studies have Identified a low hazard 
vortex wake, and these findings are presented in Section 5.2. 


5.1. Decay of a Vortex Pair 

In a calm, neutral atmosphere .- The simplest vortex wake Is 
that generated by an aircraft In cruise configuration. This wake 
Is one which can be adequately described by a simple vertex pair. 
Figure 22 shows the initial conditions for tht computation. The 
initial vorticity and turbulent kinetic energy distributions were 
taken to be Gaussian about y/s’ = ±1 . Two computations were made 
with the viscous core radius r^/s ’ = 0.4 and 0.8 . The integral 
scale parameter A/s’ was taken equal to 0.2 . The turbulent 
energy components were equally distributed with uu = vv = ww = 
q^/3 . The cross-correlations were initially taken to be zero and 
r/v = 10** . 

The computations were carried out in time and the results 
transformed to downstream distance through x = U t . Figure 23 
shows the circulation computed over the area y 0 . The nearly 
constant circulation region for the initially tight vortex 
(r^/s' = 0.4) illustrates that a time interval exists before the 
diffusion of vorticity across the wake centerline is significant. 
This time period is a measure of how long a point vortex simulation 
of the pair could be expected to be valid. 

The computation for the initially more diffuse vortex 
(r^/s' = 0.8) was carried downstream to a distance of x’/b’ = 
87rAb’^/b^CL ; at Cl = 1 and A = 7 and a rectangularly loaded 
wing , x/b 160 . For the B-747 aircraft, this is approximately 
10 km. Without ambient atmospheric turbulence, the vortex pair 
ages quite slowly. Figure 24 compares the maximum root-mean-square 
turbulence level q between the pair and a computation using a 
single isolated vortex having the same initial conditions as one 
vortex of the pair. The level of turbulence intensity in the pair 
is not significantly different from that computed for an isolated 
vortex . 

In figure 25 we show the r^/s ’ = 0.8 circulation as a 
function of downstream distance about a square of side 2e , 
centered at the centroid of vorticity in the half-plane. As might 
be expected, the decay of circulation is far more rapid than the 
half-plane value. 

Figure 26 shows the computed descent rate dh/dt of the 
vertical centroid z of the vorticity. The z increases as a 
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Figure 22. Initial vorticity and turbulent kinetic energy 

distributions used in the computation of the decay 
of a vortex pair in a quiescent atmosphere. 
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Figure 23. Comparison of the decay of circulation In a counter-rotating vortex 
pair for two different Initial vortex core radii. 




Figure 2^. Comparison of the maximum root-mean-square turbulence level In a 
^ counter-rotating vortex pair and In an equivalent isolated vortex 

u> (r_/s’ = 0.8). 
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Figure 25. Decay of circulation In a counter-rotating vortex pair computed about 

square contours centered at the centroid of vortlclty of the half-plane 

vortlcity distribution (r /s’ = 0.8). 
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function of downstream distance since a uniform upwash has been 
added to the computation (of magnitude fQ/4Trs’) to keep the 
vortex pair centered In the computational mesh. The descent rate 
decreases as vortlclty diffuses across the aircraft centerline. 

In figures 27 and 28 are shown Instantaneous streamlines at 
the beginning and end of the computation for the r^/s’ = 0.8 
computation. The y,Z coordinate system Is one In which the fluid 
at Infinity Is at rest; therefore, the pair descends downward. 

Wake dissipation by atmospheric turbulent diffusion of 
vortlclty . - There has been an attempt (ref. 31) to correlate wake 
lifetime with the turbulent dissipation rate e (our model at high 
Reynolds number gives e = 0.125q^/A). Vortices present their 
greatest hazard under stable atmospheric conditions. Under these 
conditions, scales can be smaller than the characteristic wake 
length scale — the vortex spacing. Under these problem atmospheric 
conditions, the turbulent dissipation rate does not contain suffi- 
cient Information to estimate the diffusion of trailed vortlclty In 
the atmosphere. The integral or macroscale must be known. This 
fact can be demonstrated by two computations In which e Is held 
constant. For light and light-to-moderate turbulence, we use a 
value of = 2 cm^'^^/sec as suggested In reference 31 . 

If we take A/s = 0.2 and 2.0, the corresponding turbulent 
kinetic energies become (27Tsq/r)^ = 0.025 and 0 . 116 , respectively. 
These two values can probably be achieved in a stable atmosphere 
with conventional jetliners. 

The Initial conditions are again Gaussians of vortlclty at 
y/s = ±1 of opposite sign to produce a counter-rotating pair. The 
dispersion radius Is taken to be r^/s = 0.5. The turbulent energy 
components are again equally distributed so that uu = vv = ww = 
q^/3 Initially. Integral scale parameter and turbulent energy 
components are held constant on the computational boundaries and 
r/v = 10*^ . 

Figure 29 shows the circulation decay computed about a box 
contour of side 2e centered at the centroid of the vortlclty. At 
the end of the computation (typically a wake age of 60 sec), levels 
differ by a factor of 2. More Important, the rate of decay differs 
by approximately a factor of 4. These results do well to emphasize 
the Importance of scale and caution against too simple a parameter- 
ization of the atmosphere. Obviously, the effects of wind shear 
and stratification help to make the analysis of the Interaction of 
a vortex wake with the atmosphere a most difficult one. 

An Indirect comparison can be made between the above results 
and figure 25 which gives circulation decay of a vortex pair In a 
quiescent atmosphere. The Initial spread of vortlclty Is different 
In that Vq/s = 0.8 for the pair in the quiescent atmosphere, 
while Vq/s = 0.5 for the above cases. The abscissas on figures 
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Figure 27 . Instantaneous streamlines at rt/2irs^‘^ * 0.0 and 8.0. 

The descent of the pair Is stopped by addlng^an upwash 
of r/47TS’ . Streamlines shown for rt/2Trs^ * 8.0 
have been reflected across the y = 0 axis. 
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e/s = 1.0 



Figure 29 . Decay of circulation in a counter-rotating vortex 
pair Immersed in a turbulent atmosphere with 
^ 1/3 2 
constant turbulent dissipation rate (e - 2 cm 

see figure 25. 



24 and 29 are equal, and so time or downstream positions can be 
compared directly. The slow decay of the pair in a quiescent 
atmosphere is clearly demonstrated. The dominance of the atmos- 
phere in controlling the ultimate death of a vortex wake is 
absolute even for low ambient turbulence levels. The turbulent 
atmosphere may be likened to an infinite reservoir of turbulent 
kinetic energy which can, for unlimited time, nibble away at the 
vortex wake. 

Work is underway to investigate through second-order closure 
turbulent modeling the fluid dynamics of the lower atmosphere. 

This work (partially reported in ref. 32) should give the tools 
necessary to access the rate of wake aging which results from 
interaction with a dynamic atmosphere. 


5.2. Turbulent Merging 

The result of vortex merging is the redistribution of 
vortlclty by convection and turbulent diffusion. In this subsec- 
tion, model calculations of merging are undertaken to demonstrate 
the phenomenon qualitatively and quantitatively. These results 
are then used to Identify a minimum hazard vortex wake. 

Merging computations .- An isolated vortex has streamlines 
which are circles centered on the vortex axes. The only mechanism 
which can redistribute vortlclty radially outward is diffusion, 
both laminar and turbulent. It is well known that the laminar 
aging of a vortex is an extremely slow process. Unfortunately, the 
axisyrametrlc vortex is quite stable and can only sustain low turbu- 
lence levels. In fact, turbulence is actually damped in the 
central region of the vortex (the viscous core). It is just in 
this region where vortlclty is a maximum, and the only turbulent 
redistribution of vortlclty in this region results from turbulence 
which has diffused into the core. 

That vortex cores from simply loaded wings are, in fact, 
nearly laminar can be seen in figure 5 of reference 33 where smoke 
released from a tower has been entrained into the wing tip vortex 
of a large aircraft. The air through which the aircraft flies is 
dlscernlbly turbulent (as interpreted from the dispersion of smoke) 
while the vortical fluid making up the viscous core is nearly 
laminar. 

Turbulent computations using second-order closure modeling 
also confirm this feature of Isolated vortices. Taking the equili- 
brium, nondlffusive, and high Reynolds number limits of the second- 
order closure model (ref. 2), we obtain the turbulent kinetic 
energy distribution shown in figure 30 Also shown in the corres- 
ponding distribution of vortlclty. Note that no turbulent kinetic 
energy is predicted for r/r^ £ 1.25. The centrifugal effect of 
the swirling velocity suppresses the production of turbulence in 
this region. 
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This discussion has neglected the axial velocity and, hence, 
shear as an additional source of turbulence. However, from obser- 
vations and analytic considerations (refs. 2, 3^), aircraft 
vortices do not have sufficient axial shear to provide a mechanism 
for any sizable production of turbulence. 

The presence of two or more vortices in proximity destroys 
the axial symmetry of a vortex by an induced straining field. 

Once the symmetry has been destroyed, turbulence can be produced 
which leads to the subsequent diffusion of vortlclty. As a first 
calculation of the merging phenomenon, take the merging problem 
between two like-sign, equal strength vortices. Two Gaussian 
spots of vorticlty of the form 



have been placed at y = ±d/2, z = 0 . The radius r is measured 
outward from both y = ±d/2, z = 0 . Gausslans of turbulent 
kinetic energy of the form 


q 


2 





(30) 


are also placed at y = ±d/2, z = 0 with 

2 ,2 

q^ = 0.01(r/7Td) 


The calculation is started with uu = vv = ww = q^/3 and uv = 0 . 
The integral scale parameter is initially taken equal to O.ld . 
Boundary conditions on velocity are applied via a multipole expan- 
sion of the vortlclty field (Appendix D). Turbulent quantities 
are set equal to zero at the computational boundaries. 

The results of the computations are shown in figures 31, 32, 
and 33 where instantaneous distributions of pressure, vorticlty, 
and turbulent kinetic energy are shown in intensity form. The 
counter-clockwise rotation of the pair is expected. It is inter- 
esting to note that the mean flow variables, pressure and vortl- 
city,take on a dlscernibly more axisymmetric structure at 
tr/ird^ = 4 than the turbulent kinetic energy. This is expected 
in a phenom.enon in which the redistribution of vorticlty is 
initially governed by convection. The initial time scale for 
turbulent redistribution through diffusion is of the 6'(A/q) which 
is approximately tr/7rd^ = 1 . 
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Figure 31. 


Pressure Intensity plots for the 
strength, like-signed vortices. 


merging of two equal 

p^ = - 6.0p^(r/TTd)2 . 
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Figure 32. The mean vorticlty distribution during the merg 
Ing of two llke-slgn, equal strength vortices. 
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Figure 33. Twice the turbulent kinetic energy distribution 
during the merging of two llke-slgn and equal- 
strength vortices. (qj^ird/F)* = 0.01 
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A comparison with an equivalent Isolated vortex decay compu- 
tation Illustrates the significance of the merging phenomenon In 
terms of aging the vortex flow field. We have chosen to compute 
the decay of an Isolated axisymmetrlc vortex whose initial vortl- 
city distribution Is Gaussian and whose circulation equals that of 
the pair. The core radius rc has been chosen to make the polar 
moment of the vortlclty distribution computed about the centroid 
equal to that of the merging pair. Thus, 


r 


2 

c 


r 


2 

°pair 


+ 



(31) 


The computation Is carried out with an axisymmetrlc version of the 
code (ref. 28) used to compute the merging of the like-signed pair . 
The Initial turbulence distribution is taken to be uu = vv = ww = 
q^/3 where 


2 2 
q = 2.0q^ exp 

so that the Initial total turbulent kinetic energies in both 
computations are equal. The Integral scale parameter A is taken 
to be O.ld. 

Figure 34 shows the total turbulent kinetic energy in the 
crossplane as a function of time. The level that can be supported 
by the axisymmetrlc isolated vortex is far less than initially 
Introduced, as discussed earlier and In reference 2. However, the 
breakdown of the axisymmetrlc structure about each vortex In the 
merging pair results in the production of turbulent kinetic energy 
and, hence, Reynolds stresses, which diffuse the mean vortlclty 
outward in addition to the convective’ spreading. The production 
of turbulence Is the transport mechanism which can diffuse vortl- 
clty to the wake centerline and result in the decay of circulation 
in the wake, as was shown schematically in figure 1. 

Calculations involving merging of unequal strength vortices 
have been undertaken with the second vortex 50^, 20^, and -20^ the 
strength of the first vortex and with all other initial conditions 
taken to be the same as those used in the merging calculation of 
two like-sign, equal-strength vortices. The comparison between 
the merging of equal-sign and strength vortices and merging with 
the second vortex 50% of the strength of the first can be made by 
comparing figure 32 with 35 and figure 33 with 3^. As can be 
seen, merging of equal-strength vortices involves equal redistribu- 
tion of both vortlclty distributions, while the 50 %o vortex In the 
unequal merging computation is convected about the strong vortex. 
The structure of the strong vortex Is not altered during this 
process. It Is interesting to note that while the merging of the 


'pair 


(32) 
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Figure 35 • The mean vortlcity distrihution during the merg- 
ing of two llke-slgn vortices. The vortex on the 
left is one-half the strength of the vortex on the 
right. 2Cjjr/ird* - 7.39 
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Figure 36 , Twice the turbulent kinetic energy distribution 
during the merging of two llke-slgn vortices. 

The vortex on the left Is one-half the strength 
of the vortex on the right. (qj^ird/r)^ = 0.01 
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unequal strength vortices produces a higher local level of turbu- 
lent kinetic energy, the spread of the turbulent kinetic energy Is 
greater In the llke-slgn, equal-strength computation. The results 
for the merging of the 20 % and - 20 % strength second vortices are 
not shown here, but we note that the merging redistribution Is 
qualitatively the same as observed for the unequal computation. 

The weaker vortex Is convected around the stronger vortex with 
very little change In the structure of the stronger vortex. 

A minimum hazard vortex wake .- The results given above 
may be used to determine a minimum hazard aircraft wake. The only 
assumption which Is made Is to neglect the effect of the opposite 
merging pair since the first-order effect of the vortices from the 
opposite wing Is to bring about wake descent. This, however. Is 
not true If a vortex is sufficiently close to the aircraft center- 
line so that its mate dominates. This would be the case In which 
two vortex pairs would not remain together; it has been discussed 
in Section 2.^ with the aid of figure 8. However, here we are 
Interested in merging and are looking for a configuration where 
the two vortex pairs remain together. 

If lift Is held constant in evaluating the hazard between 
various merging wakes, then 

L = pu„2y(r^ + r^) (33) 

with the positions of the vortices shown below. 



The centroid of the system Is computed to be 

af^ + sf^ 

^ ( 34 ) 
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Substituting into (33) yields 


L 

2pU 

^ oo 



l)f + 1 


(35) 


where f = . 

Vortex hazard is evaluated by NASA by measuring the rolling 
moment Induced on a following wing by the vortex wake. Although 
the calculation of aerodynamic loads on a wing using strip theory 
is of limited value even in uniform flow, it is adequate if trends 
and minimums are to be obtained. The torque Induced on a constant 
chord airfoil is 


T 



(36) 


Expressed as a rolling moment coefficient and with substitution of 
eq. (35), 


with 
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(38) 


I is the nondlmenslonal rolling moment computed in the Vortex 
Wake Code and the subscript f denotes follower aircraft. The 
decay of the maximum roj-ling moment encountered as a function of 
downstream distance for the case where the wing span of the follower 
aircraft is 20^ of the wing span of the generator is shown in 
figures 37 , 38 , and 39 . The first observation which is apparent 
from each of these figures is that the rolling moment is a minimum 
when the merging vortices are of equal strength. Selection of the 
spacing d/s is not as easy, with d/s = 0.6 and 0.8 comparing 
quite closely. However, the d/s = 0.6 spacing shows a slight 
advantage. If we refer to figure 12, the d/s beyond which equal- 
strength vortex pairs will not remain together is about 0.8, so that 
the equal strength with separation d/s = 0,6 will be a minimum 
hazard wake for a follower aircraft which has a span that is 20^ of 
the span of the generating aircraft. 

The oscillations which occur in the rolling moment as a 
function of downstream distance are expected and are associated with 
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Plffure 37. Decay of the maximum rolling moment encountered by a 
^ following aircraft of semispan 20 % of the semispan of 

the generating aircraft (d/s = 0.4) 
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Figure 38. Decay of the maximum rolling moment encountered by a 
following aircraft of semispan 20 % of the semispan of 
the generating aircraft (d/s = 0.6) 
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Figure 39. Decay of the maximum rolling moment encountered by a 
following aircraft of semispan 20 % of the semispan of 
the generating aircraft (d/s = 0.8) 




the orientation of the unmerged distributions of vortlclty relative 
to the following aircraft. As the merging process goes to comple- 
tion, the amplitude of the oscillation will decay, with the average 
value of the rolling moment (averaged between peaks and troughs) 
decaying smoothly. This can begin to be seen in figure 39 for the 
equal-strength rolling moment curve . 

If the following aircraft has a span which is 60^ of the 
span of the generator, figures ^0 and 4l show the decay of the 
rolling moment with downstream distance. Again, the equal-strength 
merging results in a minimum induced rolling moment at a given d/s. 
However, as is expected, the difference in level between various 
relative vortex strengths is not as significant as it is for the 
smaller following aircraft. The results for d/s are not shown 
here since they induce very high levels of rolling moment. There 
is little difference in rolling moment levels between the d/s = 

0.6 and 0.8 vortex separations. 

The results of the calculations given above indicate that 
the minimum hazard wake for a small follower aircraft is achieved 
by having equal-strength flap and tip vortices, with the tip 
vortices located at y = ±s and the flap vortices located at 
y = ±0.4s . This finding is a major result of the effort reported 
here. Experimental verification of it should guide future 
experimental efforts.* 

A line vortex immersed in a simple shear flow . - Moore and 
Saffman (ref. 81 have found the solution to a uniform vortlclty 
vortex in a simple shear flow. They have shown that, provided the 
simple shear is not too large compared to the vortlclty in the 
vortex, the vortex takes an elliptic cross-section. The interest- 
ing result is that when the vortlclty associated with the simple 
shear has the same sense as the vortlclty in the vortex the major 
axis of the ellipse is parallel to the streamlines at infinity. 

When the vortlclty associated with the shear is of sense opposite 
to the vortlclty in the vortex, the major axis of the ellipse is 
normal to the streamlines at infinity. 

This result was checked qualitatively, using the Vortex Wake 
Code, by placing a vortex at y = z = 0 with a Gaussian distribu- 
tion of vortlclty and superimposing a horizontal velocity equal to 
3z . The initial turbulence distribution was taken from the super- 
equilibrium limit (ref. 35) of the second-order closure model. The 
turbulent energy components are 


vv 


2uu = 2ww 



(39) 


with the cross-correlation equal to vw = -0.35A^ | dV/dz p . The 
only mean flow length scale in the problem is Vq , and the inte- 
gral scale parameter A was taken to be A = . 
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Figure 40. Decay of the maximum rolling moment encountered by a 
following aircraft of semispan 60 %. of the semispan of 
the generating aircraft (d/s ■ 0.6) 
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Results will be given for two cases Stt^Pq / F = ±0.5 . In 
figure 42 a comparison Is shown between the vortlclty of a vortex 
Immersed In a positive simple shear flow 3 < 0 and a negative 
simple shear flow 3 > 0 (the vortlclty in the shear flow is sub- 
tracted so that the level shown Is that of the Initial Gaussian of 
vortlclty). As can be seen, the results obtained by Moore and 
Saffman (ref. 8) are qualitatively substantiated. No quantitative 
comparison can be made, however, since their analysis appears to 
be tractable only when the vortlclty In the vortex Is uniform; 
numerically a uniform vortlclty vortex Is difficult to discretize. 

A comparison of the turbulence levels is shown In figure 43. 

At the end of the computation, it can be seen that the level of 
turbulence is greater In the vortex which Is Immersed in a shear 
which is in opposite sense to the vortex. This might have been 
suspected since the positive vortex is essentially Immersed in a 
negative vortlclty bath — a situation where axisymmetry is well 
known to be statically unstable. However, the issue of the produc- 
tion of turbulence by a vortex in a simple shear is more complicated 
due to the fact that sufficient shear will convect the vortex apart; 
the physical situation is far from being axisymmetric . 

Merging in aircraft wakes .- The results of three calcula- 
tions simulating the vortex interactions in the wake of aircraft 
are presented here. The first computation follows the merging 
interaction of a flap and tip vortex of like sign and equal strength. 
The tip vortices are positioned at y = ±0.95s and the flap 
vortices at y = ±0.4s. The vortices were given Gaussian distribu- 
tions of vortlclty with strength F , and = 0.2s for all 

vortices. Gaussians of turbulent kinetic energy were centered in 
each vortex with 


2 

uu = w = w = 0.033 (^) exp(-r^/r^) (40) 

and cross-correlations were initially taken to be zero. The inte- 
gral scale parameter was initially taken to be A/s = 0.2 and 
F/v = 10**. Results for the distribution of vortlclty and turbulent 
kinetic energy are shown in figures 44 and 45, respectively. 

These distributions may be compared with the results of a 
second computation where the strength of the tip vortex is reduced 
to 30^ of the strength of the flap vortex. All other initial 
conditions remain the same. This case is shown in figures 46 and 
47. 


As a third computation, a negative fuselage vortex of 47^ 
the strength of the flap vortex, positioned at y/s = ±0.2, is added 
to the initial conditions of the second computation. This vortex 
simulates the vortex trailed as a result of the lift dropoff at the 
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Figure 42. Comparison of the decay of a Gaussian of vorticlty 

Immersed In a simple shear flow having the same 

sense as the vortex with a vortex Immersed In a 

single shear flow having the opposite sense as the 

vortex (KT/-nr ~ 32 . 0 ) 
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Figure 43. Comparison of the turbulent kinetic energy between 

a vortex immersed in a simple shear of the same sense 

of the vortex with a vortex Immersed in a simple shear 

of the opposite sense of the vortex ( [q„ 2 irr./r] * - 0 . 09 ) 
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Figure ^5. Production of turbulent kinetic energy during 

the merging of the vortlclty distribution shown 
In figure 44 [q^ = 0 . 01( F^/Ptfs ) ^ ] 
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Figure 46. Distribution of trailed vorticlty in a half-plane; 

the tip vortex is 30/6 the strength of the flap 
vortex and both vortices are of the same sign 
(c_ = 50.orV27Ts*) 
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Figure 47. Production of turbulent kinetic energy during 

the merging of the vortlclty distribution shown 
In figure 46 Cq* - 0.01(rj/2¥e)*] 
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wing-fuselage Junction. The mean vorticlty and turbulent kinetic 
energy distributions for this case are shown in figures ^8 and 
respectively. 

To make some quantitative comparison between these three 
cases in terms of deintensification, the results must be scaled so 
that lift Is held constant between runs. To do this, the normal- 
ized vorticlty In figures and ^8 must be scaled up by 2.0 and 
2.3, respectively, while the nondimenslonal times In these figures 
must be scaled down by 2.0 and 2.3, respectively. The net result 
Is that the maximum level Is about one-half as great for the case 
when equal-strength vortices are merged. More significant, however, 
is how extensively the trailed vorticlty is spread by the merging 
phenomenon. The amount of shaded area In figures ^6, and 48 

is an indicator of this, with the equal-strength merging case being 
most dispersed. 

We note that In figure 49 the large level of turbulent 
kinetic energy occurs in a very small region inside the fuselage 
vortex which Is being strained by the flap vortex; It does not 
contribute significantly to the total turbulent kinetic energy in 
the crossplane. 

Engine exhaust .- It has been observed in full-scale flight 
tests that engine thrust levels do make a measurable change In the 
rate of decay of a vortex wake (ref, 36), It Is reasonable to 
suspect that the basic merging phenomenon is ndt altered by the 
presence of Jet exhausts, but that the additional turbulence pro- 
duced by the exhaust Jets results In additional turbulent diffusion 
of the trailed vorticlty. To check this conjecture, cold engine 
exhausts were added to the equal-strength llke-slgn flap and tip 
vortex merging computation described above. 

At a lift coefficient of 1.5 and assuming Cq = Ctp = 0.125 Cl, 
the axial velocity excess about each engine Is taken to be 

= 8.0 exp (- 52.0 (4l) 

Taking the aspect ratio equal to 7, the flight speed is 40.0r/27rs . 
The engine exhaust Gaussians of excess axial velocity are positioned 
at y/s = -0.75, -0.4, 0.4, and 0.75 and z/s = -0.08 . Except for 
a vertical displacement of z/s = - 0.08 , the Inboard engines 
exhaust directly into the flap vortices. Additional Gaussian dlstri 
butlons of turbulent kinetic energy are positioned at each engine 
with the distribution given by 

= 0.64 



exp - 52.0 


(42) 
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Figure 48, Distribution of trailed vortlclty in a half-plane. 
The tip vortex Is 30$ the strength of the flap 
vortex and of the same sign, A negative fuselage 
vortex Is 47$ the strength of the fuselage vortex 
Initially located at y/s = 0.2 (negative vortex 
Is not shown since black regions denote positive 
vortlclty). ” 50 . 0r^/27rs^ ) 
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The turbulent energy components are equally partitioned and the 
cross-correlations taken to be zero initially. 

Figure 50 shows a comparison of the trailed vortlclty in 
Isopleth plot form. The distributions on the left are those given 
in figure 44 with time converted to downstream distance in wing 
spans (since lift coefficient and aspect ratio are now specified). 

The smallest circles carry the largest values, since the plots 
have been made by nondlmenslonallzing the vorticity with the maximum 
value at x/b = 0 . As can be seen, the Jet engines do not hinder 
the merging of flap and tip vortices. At 12 spans downstream, the 
Jet engines give slightly more spreading of the trailed vorticity; 
hence, the lower maximum value. 

Figure 51 shows the Intensity of the turbulent kinetic energy. 
Note that the normalization differs between the run with and without 
Jet exhaust. It is likely that appreciable redistribution of the 
trailed vorticity is yet to occur owing to the large values of 
turbulent kinetic energy in the wake. This feature is overemphasi- 
zed by taking the initial integral scale parameter A/s to be 0.2 
for both runs. The Integral scale parameter in the exhaust Jets 
should be Initialized to be somewhat smaller. In any case, the 
merging phenomenon is not significantly altered by Jet exhaust and 
its accompanying turbulence. 

Jet exhaust axial velocity excess plots are shown 
in figure 52 , Note that the inboard engine exhaust, which was 
Introduced into the flap vortex, is slower to decay than the ex- 
haust from the outboard engine, demonstrating the suppression of 
the production of turbulence by swirl. Engine exhausts are, of 
course, not the same temperature as the surrounding fluid. The 
density difference between the exhaust fluid and the atmosphere 
gives rise through gravity to a buoyancy body force on the exhaust 
which becomes more significant with wake age. To demonstrate the 
buoyancy effect, a sample computation was undertaken of four hot 
Jets exhausting into an atmosphere. The engines were positioned 
as described above. No vortex wake was Initially specified so that 
the generation of vorticity and turbulence by buoyancy could be 
more easily discerned. A temperature profile was assumed at the 
exit of each Jet which would give rise to a strong buoyancy effect 
(to check the code): 


AT 

T 


= 0.5 exp ( - 52 . 0 — 5 - 


( 43 ) 


The heat added by this distribution is rather larger than would be 
added by an existing Jet engine. Keeping total thrust equal to the 
above computation requires 


AUtts = i 4 .o exp (- 52.0 


( 44 ) 
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Figure 52. laopleth plots of the Jet exhaust axial velocity 
excess. = 8.0r/2irs . 
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The turbulent kinetic energy distribution about each engine v/as 
given by eq . (^2 ) . 

The resulting decay of the exhaust velocity is shovm in 
figure 53* Without the presence of the vortex swirl velocity and 
distribution of axial velocity excess, the decay is quite different 
from that shown in figure 52, In figures 5^ and 55 are shown the 
distributions of temperature and turbulent kinetic energy. Mote 
that the normalizing value for each distribution in figure 55 is 
changed to allow intensity resolution over the large change in 
turbulent kinetic energy. The vertical upward motions of the 
distribution are evident in figures 53? 5^, and 55* 

The axial or streamwlse vorticity produced by buoyancy is 
not easily presented in intensity plot form since both positive and 
negative vorticity are generated (the total circulation is, however, 
zero). Instead, in figure 56 the maximum and minimum values of 
vorticity are shol^m as a function of downstream distances. 

V/hile the area integral of the total vorticity generated by 
buoyancy is zero, the y moment which is associated with the 
vertical impulse given to the fluid by the body force is not. Since 
this negative impulse in a neutral atmosphere is ever increasing at 
some time after aircraft passage, the total impulse in the wake is 
zero. Estimates of this time have been difficult to make vjith any 
certainty, but minutes rather than hours appear to be the correct 
order of magnitude estimate. V/hile admittedly the amount of engine 
heat used in the above calculation is larger than that discharged 
by a conventional by-pass ratio engine, it is easy to suspect that 
buoyant effects as well as the turbulence produced by the kinetic 
energy of the jet will explain the sensitivity of vjake decay to 
thrust level settings. 
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Intensity plots of the hot jet exhaust axial 
velocity excess AU^ = l4.or^/2xs 
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Figure 55* Intensity plots of the turbulent kinetic energy 
corresponding to the axial velocity excess shown 
in figure 53 
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6. CONCLUSIONS AND RECOMMENDATIONS 


Several conelusions from this study can be made with regard 
to predicting aircraft vortex wake strength and structure. 

1. The number, strengths, and initial positions of vortices 
trailed from an aircraft wing are determined by the detailed lift 
distribution on the wing. Prior to this study, only Invlscld models 
were available to predict the swirl velocity distributions in each 
vortex and the Interactions between vortices. The vortex-vortex 
interaction, called merging, which has been shown in the laboratory, 
in flight tests, and through analytic techniques to lead to rapid 
aging of a wake involves the turbulent diffusion of trailed vortl- 
city. A transport model is now available to predict viscous wake 
interactions . 

2. A review of invlscld theories to predict wake structure 
suggests that the Betz technique is only of limited usefulness in 
multiple-vortex wakes where strong interactions occur between 
vortices . 


3. Point vortex computations using large numbers of line 
vortices to represent the initial sheet of trailed vorticlty and 
the subsequent motion of these vortices provide useful information 
during the initial stages to wake development. This technique, 
however, excludes the mechanism of viscous diffusion and, therefore, 
conserves circulation and fluid kinetic energy in the wake half- 
plane. Extending the model to Include the effects of engine exhaust, 
atmospheric turbulence, and buoyancy is not straightforward. 

4. Point vortex computations of the rollup and interactions 
of the vortices making up the wake of a B-747 aircraft have verified 
the trends in terms of Induced rolling moments on a following air- 
craft that have been observed in the laboratory and in flight tests. 

5. A study using point vortices has determined that the 
merging of flap and wing tip vortices can be sensitive to small 
changes in fuselage vortex strength and position. 

6. A study using point vortices has determined that the 
merging of flap and wing tip vortices can be sensitive to the proxi- 
mity of the ground. 

7. Calculations using the vortex wake transport code 
developed during this study have shown that a minimum hazard wake 
is trailed from a generating aircraft when the aircraft configura- 
tion is such that a flap and wing tip vortex of the same sign and 
strength are shed from each wing. To achieve this optimum, the flap 
vortex is located outboard approximately 40^ of the distance to the 
tip vortex. The vortices merge, during which convection and diffu- 
sion distribute the trailed vorticlty widely over the wake. 
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8. Calculations using the vortex wake transport code have 
shown that engine exhausts do not appreciably alter the merging 
phenomenon. Results suggest that the turbulence added to the wake 
by engine exhausts help further to diffuse the trailed vortlclty 
prior to, during, and after vortex merging has occurred. 

9. The effect of ambient atmospheric turbulence on the 
aging of an aircraft wake was Investigated at constant turbulent 
dissipation rate. It was shown that under stable atmospheric condl 
tlons, when atmospheric macroscale may be less than or equal to the 
vortex spacing, misleading results may be obtained If vortex aging 
Is correlated only with the turbulent dissipation rate. This 
result cautions against using one parameter to characterize the 
ability of the atmosphere to dissipate aircraft wake vortices. 

The results obtained from this study warrant the following 
recommendations . 

It Is becoming Increasingly apparent that the key to being 
able to explain why small changes In the configuration of an air- 
craft can lead to a substantially stronger or weaker wake lies, to 
a large extent, in knowing what has been done to the wing lift 
distribution. Analytical techniques either do not exist or appear 
to have dlfflculy In predicting the lift distribution on a wing 
with landing gear, flaps, and spoilers deployed. It is recommended 
that accurate load distributions be obtained experimentally for the 
B-7^7 aircraft, particularly for configurations whose wakes have 
been probed or measured. The load distribution data will permit 
explanation of existing data showing configuration sensitivity and 
will permit minimization techniques to be applied to other aircraft 

Results of studies presented here suggest that the vortex 
merging interactions in the vicinity of the ground should be the 
subject of experimental and analytical studies. A minimum hazard 
configuration cannot be strongly altitude- or aircraft attitude- 
sensitive . 

Methods to assess the ability of the atmosphere to dissipate 
wakes are now available. In light of the fact that measurement of 
wake decay In a real atmosphere Is an enormously difficult task, 
these methods should now be exploited. 
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APPENDIX A. 


ANALYTICAL PREDICTION OF MERGING BOUNDARIES 


Moore and Saffman (ref. 8) have given an analysis of a vortex 
of elliptic cross-section in which the vortlclty is uniform and 
have determined far field flows in which steady, stable solutions 
are possible. The results of this analysis are extended to deter- 
mine the proximity d of a vortex of strength T 2 to a vortex of 
initial diameter Uq and strength Fi (uniform vortlclty) which 
will not yield a steady solution of elliptical cross-section. The 
two vortices are shown in the sketch. 



A solid body rotation flow field is added to cancel the rotation. 
The flow field is then steady. The streamfunction for the solid 
body flow field is 


= I" f2[(x - x)^ + y^] (A.l) 

2 — 

where ^2 = (P-, + r2)/27Td and x is the centroid and is given by 
X = -d/(i + ri/r2). 

The flow field Imposed on the vortex at x = y - 0 is given 
by 

^ (x^ - y^) + ~ fi (x^ + y^) + e^(xVd^) (A. 2) 

Substituting into Moore and Saffman’ s equation (2.l6C) with 
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r 21 

e = — ^ , y - y and w = — \ - 2Q 


2Trd 


^ TTd ^ 
c 


yields 


0^(s + 2) + 0^(-2st^ + s + 2) + 0(2st^ - s) - s = 0 (A. 3) 


where s = ^ “ 2d/d(j . Real positive roots with 

0 = a/b > 1 have the major axis of the ellipse on the x-axis (a/b 
is the major to minor axis ratio). Those roots lying 0 < 0 < 1 
have the major axis along the y-axis. 

The solutions obtained by solving eq. (A. 3) must be checked 
for stability using Moore and Saffman’s equation (3.17) • The 
curves along which steady, stable solutions cease to exist in terms 
of the parameters ^/^c shown in figure 3- 
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APPENDIX B. 


COMPUTATION OP A MEASURE OP THE ROLLING MOMENT INDUCED 
ON A POLLOWING AIRCRAFT IN THE WAKE OF A 
Discrete Vortex Wake Model 


First, the wing load distribution for the aircraft configura- 
tion being studied is needed. A standard approach to predict 
this distribution analytically is vortex lattice theory. The wing 
is modeled by a lattice of horseshoe vortices whose strengths are 
such that the flow normal to the mean camber line of the wing at 
specified points is zero; i.e., no flow through the wing is 
allowed. The lift associated with each vortex is computed and 
summed appropriately to determine the spanwlse wing lift distribu- 
tion, total lift, and pitching moment. From the latter two 
quantities plus the c.g. position, the tall lift required to produce 
equilibrium flight (no pitch acceleration) is calculated. 

Then, a desired number of constant strength line vortices, 

2N , is positioned along the span so as to correspond to the pre- 
viously calculated wing load distribution, where N is the number 
of vortices for each half of the wing. (One of the N pairs of 
vortices may be of lesser strength since the spanwlse load is not 
necessarily divisible by an Integer number of Incremental loads.) 

An elliptic load distribution is assumed for the horizontal tall, 
and the resulting single vortex pair is located laterally at 

TT ^T 

y = ± jj- where b<p is the tall span. The vertical position of 

this vortex pair is determined by aircraft geometry and angle of 
attack. The total number of vortex pairs comprising the wake is, 
thus, N + 1 . The trajectories of these N + 1 vortex pairs as 
they convectively interact with each other are calculated in time 
or, equivalently, in downstream distance (see fig. Bl). 

With the location of the N + 1 pairs of vortices given in 
the cross-flow plane as a function of downstream distance x , 
there is yet the problem of determining how hazardous the wake is 
to an encountering aircraft . This could be approached by using 
strip theory to determine the roll torque as a function of aircraft 
position and orientation relative to the wake. Unfortunately, the 
unrealistically high velocities near the center of each point vortex 
make this approach singular whenever a vortex lies on the follower’s 
wlngtlp. This could be circumvented by postulating a vortex core 
region, but then the results would be sensitive to the core size 
chosen. A much simpler approach is to assume that the circulation 
computed in an area characteristic of the dimension of the encount- 
ering aircraft is some measure of hazard. One could then simply 
survey the wake with a square or circular window with dimension 
equal to the follower aircraft span and sum the circulation in the 
window. This approach, while simple, has unrealistic step changes 
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centerline 


Figure Bl. Geometry for the calculation of the downstream 
location of vortex centroids. 


as vortex centers cross the boundary of the window. Therefore, a 
more physically motivated criterion Is selected which has no mathe- 
matical problem with the singularities at the center of each point 
vortex. The quantity chosen as the basis for this criterion is the 
angular momentum within a circle of diameter equal to the wingspan 
of the follower aircraft. For it is this angular momentum in the 
fluid surrounding the encountering aircraft that can be transferred 
to it as a roll torque or rolling moment. Thus, an aircraft roll 
rate is directly related to the angular momentum of the surrounding 
fluid. 


The mathematical details of these procedures, as depicted in 
figure 4, to compute the roll hazard in the wake of a generating 
aircraft are given below. While the method is general, the details 
are given specifically for the B-74? aircraft. 

Vort ex lattice method for wing load distribution .- The first 
step in the analysis is to compute wing geometry of the B-747 as a 
function of trailing edge flap and leading edge slat positions. In 
particular, the local angles of attack of the mean camber line, 
must be specified at a series of points on a grid covering the wing. 

The grid locations and angles of attack are functions of the 
basic wing planform, twist, and airfoil sections, and the positions 
of the leading edge slats, ^^leF » Inboard and midspan 

trailing edge flaps, and . The wing planform of the 

B -747 with and without surfaces extended is shown in figure B2, The 
mean camber angle at any spanwise location is given by linearly 
interpolating between the angles given as a function of x/c in 
Table C .1 (Appendix C ) where x/c is the chordwlse station from 
the leading edge nondimenslonallzed by the chord length at that 
spanwise location. In addition, the outer panel of the wing is 
twisted with a linear spanwise variation reaching - 3-1/2° at the 
tip . 


^“twist ^ -^Ky’ - .71) 


y’ 1 .71 


y’ 


= y 

b/2 


(B.l) 


The major changes to local angle of attack and planform occur 
when the flaps and slats are deployed. For the purpose of this 
study, both devices are assumed to be uncambered surfaces with 
sealed slots. The surfaces extend from and pivot about a fixed 
hinge line, thus positioning the endpoint of the surface a horizon- 
tal distance Ax and vertical distance Az from the pivot point 
as Illustrated in figure B3. Also, in keeping with the accuracy of 
this analysis, the three leading edge slats are approximated as a 
single surface with no spanwise variation. 
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Figure B2. B-7^7 wing planform with and without flaps and slats 
extended (data supplied by The Boeing Company) 
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The relations between flap or slat angle setting and these 
two coordinates (Ax, Az) are given In tabular form in Appendix C. 
Subroutines have been developed to Interpolate among these data 
and thereby to specify planform and camber angle changes necessary 
to predict the wing load distribution. 

It then remains to divide the wing Into panels, each of 
which Is represented by a horseshoe vortex. Each panel produces a 
lift at 1/4 of the panel chord measured from the leading edge. In 
response to the downwash velocity computed at 3/4 of the panel 
chord. A grid containing l80 panels used for this study Is shown 
In figure B4. The number of equal length panels In the chordwlse 
direction is fixed at each spanwlse location chosen. The actual 
length varies with flap and slat deflection. 

Thus, for a given set of slat-flap deflections, a mesh of 
values of local angle of attack is computed and used as Input to a 
modified version of Margason's vortex lattice program (ref. 11), 

Rather than prescribing a wing lift coefficient, the program 
has been modified to accept total lift coefficient and to compute 
the tall and wing lift contributions. This is done using the 
equilibrium lift and pitching moment equations as shown below. 



Lift: 




(B.2) 


Moment about 
c .g. : 

where 


X, 


‘Cp^ ^cpm 

-‘-'T 


o 
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= lift coefficient 


o 

S = 
c = 


X 


= moment coefficient at = 0 
area 
chord 

location of center of pressure 
wing 
T = tall 


cp 

W 


(B.3) 



Total panels = l80 at 20 spanwise stations 

(The chordwise dimension of the panels at any 
spanwise station on the flaps is a function of 
flap position ) 


Figure b 4 . Subdivision of wing into panels for calculation of 
wing load distribution by vortex lattice theory 
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Combining eqs. (B.2) and (B.3) to eliminate 0^. yields 


CtX - 

L cp^ W 






cp 


- X 


(B.4) 


T 


cp 


W 


For this study , the total lift coefficient is Cl = 1.2 . 

Using a constant lift curve slope, the aircraft angle of 
attack a is calculated, and the wing load distribution is obtained 
as 


where 




= f(y') 


c^ = sectional lift coefficient 
c = local chord 


C-r = total lift coefficient 
c^ = reference wing chord 

In addition, the tail load magnitude is calculated from eq. (B.2). 
c^c 

These data, — ^ (y’)s hm 5 and a are then used to generate a 
ClCw = 

point vortex model of the flow field. 


"Discretization" of the wing load distribution and computa- 
tion of the point vortex aircraft wake .- Given the wing load 
/coC \ 

distribution | at locations y. , z. and the tall load 

\CLCw/i 

at ytp , Zm , which have been calculated for a particular flap-slat 
configuration, it is necessary to model the aircraft wake flow 
field using a finite number N of vortex pairs for the wing and 1 
for the tall. For the computations undertaken here, 30 £ N £ 99, 
with approximately 40 used most of the time. The discretization of 
the trailed vortex sheet calculates the locations of vortices of 
constant strength as follows. 


Ac ^c 

An initial Incremental load, — — = Ail^ , is chosen and 

Cl=w 

starting from the wing tip, the file is checked until a value 

CpC 

^ Is found. The exact location at which — — = 

1 



3 




o 
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Is found by interpolation. The search then resumes in the inboard 
direction until an Increase or decrease in load from nA^Q (in this 
case n = 1) of AJIq is found. This process is repeated across 
the semispan of the wing until y* = 0 is reached. At this point, 

/ \ 

a vortex of reduced strength corresponding to ( — ^ - nAil 

\^L°W/y» = 0 

is computed at the spanwise location y’(nA£g)/2 . The number of 

vortices computed by this process is compared with the desired 

number and, if , the initial Incremental load 

’ computed desired ’ 

is Increased by 10%, A£i = 1.1 A£q, and the process is repeated 
until ^computed — ^desired ’ 


In addition, the vortex strength associated with the tail 
load is computed as 


r 


T 






(B.5) 


The tall loading is assumed elliptical so that 


y 


T 


y ^ 

W2 " ~W 


(B.6) 


The vertical location of this vortex relative to the wing vortex 
plane is a function of angle of attack as shown below. 



The wing vortex strengths and initial positions ^i a yi > 
z-^ , where = 0 and, similarly, for the single vortex tail pair 

, yip , hip , are stored, and the calculation of their motion is 
undertaken. 
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The velocity of a point vortex in the wake field is obtained 
by summing the ■ velocity contributions made by all the other 
vortices at the center of the vortex In question. This is 


y = U 

V rn 


dy 


N+l 


m _ 


m dx 


2tt 


1=1 

l?^m 


- ( z -z . ) 
m 1 


z -z . 

m 1 




(B.8a) 


dz. 


z = U 


m 


m dx 


N+l „ 
2ir 

1 = 1 
Ij^m 




V=^i 


(z ~z.)^ + (y -y.)^ 
^ m 1 '^m 1 


(z -z . ) ^ + (y +y . ) ^ 
m 1 *^m ‘^1 


m 


Wy 


(B.8b) 


m 


where positive circulation T following the right-hand rule is 
represented by a vector in the +x direction. Nondlmensionalizlng 

X by C^/bnA where A = wing aspect ratio 

y by b/2 

z by b/2 

and substituting for 


r 

1 


(^) 


Uc..Ct 
W L 


(B.9) 


the trajectory equations become: 


dy ' 
dx’ 


N+l 



- ( z ’ -z ! ) 
m 1 




+ 


( y ’ ~yi ) 


z ’ -z 
m 


(z ’-z ! )^ + (y ’+y ! )^ 
m 1 ‘^m 1 


(B.lOa) 


dz ’ 
m 

dx’ 


N+l 



y * ~y • 

•^m 1 


(z;- 4)2 + (y;-yp2 




( z’ -z’. ) ^ + (y’ +y’. ) ^ 
^ m 1 ^ ^ m If 
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(B.lOb) 


m 


2 y ’ 


These 2N + 2 first-order differential equations are integra- 
ted using a modified Runge-Kutta procedure, and the locations , 

of each vortex saved at specified values of x’ 

In addition to specifying the numerical procedure used, a 
word is in order about the Integration step size versus accuracy. 
Since the system of equations is large, it is particularly neces- 
sary to use a small enough step size to track the motion accurately 
but not small enough to incur significant roundoff errors . 
Fortunately, the physics of the process requires that the Impulse 
and kinetic energy associated with the relative position of the 
vortices (Klrchhoff-Routh path function (ref. 2)) in the wake model 
remain constant. Therefore, these quantities can be monitored to 
determine if, and check when, the integration process is signifi- 
cantly in error. In general, however, the Integration process 
could be carried out to a nondimensional downstream distance of 
unity (about 1.5 km for the B-747) with only a few percent change 
in the invariants. 

This completes the procedure for defining the aircraft wake 
flow field in terms of trajectories of invlscld point vortices. 

The remaining step is to evaluate how hazardous the flow field is 
for smaller following aircraft. 
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APPENDIX C. 


TABULAR DATA FOR BASIC WING CAMBER AND FOR PLANFORM AND CAMBER 

CHANGES DUE TO FLAP EXTENSION 


The manufacturer supplied mean camber data for the wing with 
no flap extensions, as given below, as a function of nondimensional 
chordwlse position x/c . The characteristic camber is the same 
for all spanwise locations. 


x/c 

°^local 

(rad) 

0 

-.2472 

.025 

-.2034 

.050 

-.1628 

.075 

-.1259 

.10 

-.0926 

.15 

-.0381 

.20 

.0001 

.25 

.02168 

1.00 

.0270 


TABLE C.l 



Since both leading and trailing edge flaps are not actually flat, 
the angle variable 6 does not geometrically describe their angu- 
lar position. It is rather a nominal value used to specify the 
positions corresponding to linear extensions Ax , Az as ^listed 
in the following tables and illustrated in figure B 2 . Points 
falling between data points are obtained by linear interpolation. 
For the midspan flap, it is also necessary to interpolate in the 
y direction, using the data given at 2 y/b = O.63 and O.50. 
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TABLE C.2 


Leading; Edp;e Flaps 


'^LEF 

0 

60 


Ax (cm) 
0 

46.48 


Az (cm) 
0 

66.04 


Trailing; Edg;e Flaps 


Inboard Flap 





6^ (deg) 

Ax (cm) 

Az (cm) 



0 

0 

0 



2.5 

108.71 

14.98 



5.0 

129.28 

30.22 



H 

CO 

o 

154. 9^ 

81.13 



25.5 

161.54 

208.21 



28.67 

184.91 

234.18 



Midspan Flap 






^^.63 


''^56 

''^56 

0 

0 

0 

0 

0 

2.5 

61.36 

10.46 

61.26 

11.68 

5.0 

93.26 

20.62 

93.37 

23.06 

8.5 

99.87 

31.29 

99.77 

35.30 

14.75 

109.52 

55.77 

108.71 

62.99 

19.25 

120.09 

95.19 

110.94 

107.44 

23.0 

124.15 

137.36 

111.76 

155.19 
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APPENDIX D. 


NUMERICAL BOUNDARY CONDITIONS FOR A FLUID OF INFINITE EXTENT 


A new technique Is devised to specify computational boundary 
conditions when the vortlclty Is not extensively distributed 
throughout an Infinite two-dimensional fluid. Consider the 
vortlclty distribution sketched below at time t . 



The fluid Is Incompressible and the motion two-dimensional. At 
time t the stream function Is computed from a Poisson 

equation 

= -C(y,z,t) ^^^2) 

The solution to (D.l) is 


^(y,z,t) 





>t) 


log ((y-y’)^ + (z-z’)^^dA' 

(D.2) 


where the logarithm Is the appropriate Green’s function. 
The V and w velocities are obtained from 


W 


ii 

3y 


V = + 


3i|; 

3z 


(D.3) 
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Now, in principle, the integral in eq. (D.2) could be carried out 
numerically and the velocities could be obtained at the computa- 
tional grid boundaries at time t . If the computational grid 
contained NM points, on the order of 2(N + M) Integrations over 
the computational domain would need to be carried out. However, 
the computer time required would be prohibitive. Therefore, an 
alternate approach is necessary. 

If we fix our numerical computational boundaries so that they 
are well away from and surround the region containing vortlcity, 
eq. (D.2) may be approximated by expanding for y >> y^ and 
z >> z' . The result is 


where 


and 


00 00 


i^(y,z,t) = - 





mn mn 


‘mn 


/ oo 


m=0 n=0 


m 


C( y * , z ’ , t ) dA 


°oo = 


G 


01 


G 


10 


G 


02 


G 


11 


G 


20 



(D.4) 

(D.5) 


The Integers n and m are sufficiently large to give the desired 
resolution of the far field. Our work suggests that n and m 
need be no larger than 2, The velocity boundary conditions computed 
from the vorticity at time t are used at time t + At , giving 
first-order accuracy. One iteration could be used to raise the 
boundary condition to second-order accuracy, but this has proven to 
be unnecessary in the Vortex Wake Code a^t this time. 


Some additional economy can be gained by suitably choosing 
the origin of the y,z system to coincide with the centroid of the 
vortlcity distribution. The centroid is defined by 
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y ^01^^00 

z = 

which, obviously, is meaningful only if Iqo ^ ^ • When this is 
done, the velocity decays as r where r is measured from the 
centroid. The next term decays as r . 

When the wake of an aircraft is computed, a half-plane calcu- 
lation can be used. In this case, an image system is used to 
obtain the induced field from the reflected plane as shown below. 



The details of the superposition of velocities are straightforward. 
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